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Abstract 

In this paper, we investigate the role of a high-frequency magnetic field in the resonant behavior displayed by a 
spin- 1/2 particle under the influence of a rotating magnetic field. We propose two alternative methods for analyzing 
the system dynamics, namely, the averaging method and the multiple scale method. The analytical results achieved 
by applying these two methods are compared with those obtained from the numerical solution of the Schrodinger 
equation. This comparison leads to the conclusion that the multiple scale method provides a better understanding of 
the system dynamics than the averaging method. In particular, the averaging method predicts the complete destruction 
of the resonant behavior by an appropriate choice of the parameter values of the high-frequency magnetic field. This 
conclusion is disproved both by the numerical results, and also by the results obtained from the multiple scale method. 

Keywords: Coherent destruction of tunneling, resonant behavior, spin-1 /2 particle, rotating magnetic field, 
high-frequency fields, two-state systems 



1. Introduction 

The dynamics of quantum systems under the influ- 
ence of time-periodic driving fields is a long-standing 
problem in Quantum Mechanics d 0, S H H]. An 
extensive review on this topic, covering both non- 
dissipative and dissipative systems, can be found in 
Ref. ||6[]. For a more concise review with an accurate 
description of the Floquet-theoretical method, it is rec- 
ommended to consult chapter 5 in Ref. IM ■ 

In the past two decades, the study of driven quantum 
systems has sparked renewed interest after the discovery 
of the phenomenon of coherent destruction of tunneling 
||8[]. In the original work ||8|], the authors considered the 
problem of a particle moving in a double-well potential 
and perturbed by a monochromatic driving field. The 
analysis of the Floquet spectrum of the system allowed 
them to show that, by an appropriate choice of the pa- 
rameter values of the driving field, the tunneling of the 
particle through the potential barrier can be brought to 
a standstill. Among the possible applications of this 
phenomenon, one can mention the laser-induced trap- 
ping of an electron in a quantum- well structure Jgl llOll . 
and the control of proton or electron transfer reactions 
illi. 12, 131 . Recently, the phenomenon of coherent de- 
struction of tunneling has been experimentally verified 



in an optical double well system f 14], in a Bose-Einstein 
condensate in a strongly driven optical lattice 1 15] and 
in single-particle tunneling in strongly driven double- 
well potentials llal . 

In order to clarify the mechanism of coherent destruc- 
tion of tunneling reported in Ref. [8J, it has turned out 
to be useful to consider a simplified two-state model 
lISHHIlSIIIllIlllllMl- In this model, all 
the spatial information contained in the Floquet modes 
is neglected, and only the influence of the lowest 
quasienergy doublet is taken into account. The study 
of this approximate two-state model makes it possible 
to determine simple conditions for the appearance of 
the phenomenon of coherent destruction of tunneling. 
For instance, in the high-frequency limit, an averag- 
ing procedure shows that this phenomenon occurs when 
the ratio of the amplitude of the monochromatic driving 
field to its frequency is a zero of the zeroth-order Bessel 
function of the first kind 112 ill . It is worthwhile to point 
out that the phenomenon of coherent destruction of tun- 
neling in a two-state model resembles that of dynamic 
localization observed in a driven, infinite tight-binding 
system |25j. A discussion of their similarities and dif- 
ferences can be found in Sec. 4 of Ref. |6J. Recently, 
the internal relationship between these two phenomena 
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has been clarified in Ref . ll26ll . 

In the last decade, several works have been published 
which describe alternative approaches, beyond the av- 
eraging procedure, for the analytical study of a driven 
two-state quantum system |27, 28, 29l l30ll3U 13211 . In 



particular, in Refs. 1.27. .28. .29.1 the authors propose 
a convergent strong- coupling expansion free of secu- 
lar terms. In Refs. [30l blf the method is based on 
dual Dyson series and renormalization techniques. Fi- 
nally, in Ref. ll32ll the author applies a multiple scale 
expansion in the inverse of the driving frequency for 
the calculation of the Floquet states and the quasiener- 
gies. One of the main conclusions of these works is that, 
for high but finite values of the driving frequency, the 
coherent destruction of tunneling in a two-state quan- 
tum system is only temporary. More precisely, in the 
high-frequency limit, coherent destruction of tunnel- 
ing occurs only for time values much smaller than a 
characteristic time scale proportional to the driving fre- 
quency squared. This is a consequence of the fact that 
the quasienergy degeneracy predicted by the averag- 
ing procedure when the ratio of the amplitude of the 
monochromatic driving field to its frequency is a zero 
of the zeroth-order Bessel function of the first kind is 
exact only in the mathematical limit of infinite driving 
frequency. For high but finite values of the driving fre- 
quency, there is a quasienergy splitting of the order of 
the inverse of the driving frequency squared (see, for 
instance, Ref. [32']). 

In the present paper, we take up again the problem 
of a spin- 1/2 particle in the presence of a rotating mag- 
netic field, originally studied by Rab i 1211 . This same 
model has been considered in Ref. 13311 to study the 
suppression of quantum coherence induced by circu- 
larly polarized driving fields. As is well known, under 
the influence of a rotating magnetic field, the dynamics 
of a spin- 1/2 particle displays a resonant behavior (see, 
for instance. Sec. 11.2 of Ref. i34ll or the beginning of 
Sec.|4]of the present paper). In this work, we are inter- 
ested in analyzing how this resonant behavior is affected 
by the application of an additional high-frequency mag- 
netic field along the rotation axis. As will be shown in 
this paper, there is a remarkable similarity between the 
problem considered here and the problem of coherent 
destruction of tunneling in a two-state quantum system 
driven by a monochromatic field. 

The paper is organized as follows. In the following 
section, we introduce the model and apply an averaging 
method to work out an analytical expression for the time 
evolution operator In Sec.|3] we use multiple scale tech- 
niques 135, 361 to obtain an improved expression for the 
time evolution operator which is valid in a wider range 



of time scales than the one obtained by the previous av- 
eraging procedure. In Sec.|4] we apply the analytical re- 
sults in the previous sections to the study of the effects 
of the high-frequency magnetic field on the dynamics of 
the system. In particular, we examine in detail its effects 
on the above mentioned resonant behavior. In that sec- 
tion, we also compare our analytical results with those 
obtained from the numerical solution of the Schrodinger 
equation. Finally, in Sec.|5] we present conclusions for 
the main findings of our work. 

2. Description of the model. The averaging method 

The effect of a uniform, gyrating magnetic field of the 
form 



B(f) - Bj_ cos(ajt)ex + sin((yf)ey + Byez 



(1) 



on the dynamics of a spin- 1/2 particle has been well- 
known since the classic paper by Rabi [2J. In the above 
expression ex, ey and Ci are the 3-D unit vectors along 
the X, y and z axes, respectively, the constants By and Bj^ 
are the components of the magnetic field parallel and 
perpendicular to the z axis, respectively, and w is the 
angular frequency of the magnetic field's rotation [see 
Fig. O]. 

The time evolution of a state vector I*?, t) is governed 
by the Schrodinger equation 



d\^ t) 



(2) 



where y is the gyromagnetic ratio, S - h{&^ex + cTyCy + 
&,ez)/2 is the 3-D spin operator (with &x, o-y and &, 
being the standard Pauli matrices), and the centered dot 
denotes the usual scalar product of two 3-D vectors. 

One of the main features of this system is the fact 
that, under the influence of the rotating field in Eq. ([T]i, 
the spin dynamics displays a resonant behavior This 
resonance behavior manifests itself in several ways. For 
instance, as discussed in Sec. |4] if one takes as initial 
condition I*!*, 0) an eigenstate of &- and assumes that 
Bj_ + 0, the expected value (o-j, t) - (^, t\&~^, t) is an 
oscillating function of time with maximum amplitude at 
an angular frequency w = -yB\\. 

The aim of this paper is to study how this resonant be- 
havior is affected by the presence of an additional high- 
frequency magnetic field of the form 



BHF(f) = BhF cos (i^Hpf + f) Cz, 



(3) 



where Bhf is the amplitude of the high-frequency field, 
Qhf is its frequency, and (p e [0, 27r) is an arbitrary 
phase shift between Bhf(0 and B(f). In this case, the 
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Figure 1 : Sketch of the gyrating magnetic field in Eq. (T). 

equation of motion of the system is given by Eq. ^ 
with B(f) replaced by B(f) + Bhf(0- By high-frequency 
we mean that Qhf is much higher than the rest of the 
characteristic frequencies of the problem in the absence 
of Bhf(0- Thus, if we introduce the characteristic fre- 
quencies ct)|| = yB|| and a>j_ = yB^, it is assumed 
that a>^, (jj||, CO <& Qhf- In principle, the frequency 
whf - jBhf may be of the same order of magnitude 
as Qhf- In what follows, lj will be taken as the unit of 
frequency (i.e., a> - 1), so that cot = t, co±/a> - coj_, 
a>\\/a> — W||, whf/w = whf, and Hu^/o) — Qhf are all 
dimensionless quantities. 

In this section and the following one, we are inter- 
ested in obtaining analytical expressions for the time 
evolution operator from time f = to time t. This oper- 
ator, which will be denoted by U{t), fulfills the differen- 
tial equation 



with the initial condition 

t/(0) = /, 
where 



(4) 



(5) 



H(t) = — -i^ [cos(Oo-x + sin(Oo-v] 

--[oJ\\+ whf cos (Qhf? + ip)] o-. (6) 

is the Hamiltonian operator in dimensionless units, and 
/ the identity operator. 



The differential equation (HJi can be transformed to a 
more convenient form for the calculations by perform- 
ing the following transformation 



G{t) = e' 



At-eumji 



U{t), 



with 



d(t) - rsm(D.uFt + <p). 



(7) 



(8) 



and r - whf/^hf- Then, it is straightforward to see 
that the transformed time evolution operator G(t) can be 
obtained by solving the differential equation 



dG(t) . 
ot 
with the initial condition 



G(0) = e 



-i0(O)&:/2 



(9) 



(10) 



and the transformed Hamiltonian 



/z(QhfO = 



-^{cos [0(0] (r, + sin [0(f)] <Ty} 



1-HWii 



(11) 



For the multiple scale method developed in the next 
section Eq. (|9]l offers two advantages overEq. (HJl. First, 
in contrast to the original Hamiltonian H(t), the trans- 
formed Hamiltonian /!(QhfO does not contain terms 
proportional to whf which, as mentioned above, may be 
of the same order of magnitude as the high-frequency 
Qhf s> 1 . Second, as we have made explicit by using 
the notation h{Q.iiFt), the transformed Hamiltonian de- 
pends on time only through Qhf? which, as we will see 
in Sec. [3l is the rapid time scale to in the multiple scale 
analysis. 

In the absence of the high-frequency field (i.e., for 
r - 0), Eq. (|9]l can be solved exactly since the trans- 
formed Hamiltonian in Eq. ( fTTT i becomes time indepen- 
dent. This is not the case for r ^ and, consequently, 
one has to resort to approximate techniques in order to 
obtain an analytical solution of Eq. ^. In this section, 
we will apply an averaging method to solve Eq. (|9j. 
From Eqs. (|9]l and ( fTTT i it follows that the time deriva- 
tive of G{t) is at most of order 1 in the high-frequency 
limit Qhf » 1 ■ Therefore, a large number of oscilla- 
tions of the function 0(t) appearing in the definition of 
/i(QhfO takes place before a significant change in G(t) 
occurs. As a consequence, for Qhf » 1, the time de- 
pendent Hamiltonian /!(QhfO appearing in Eq. (|9]l can 
be approximated by its average over one high-frequency 
period. The result is 



I—:: — = «effG(0, 
ot 



(12) 



where hes is the effective Hamiltonian 



leg 



In 



■27T 



dtQ h(tQ) 



= -J [co±JQ(r)&j + (1 + W||)o-e 



(13) 



with Jo(r) being the zeroth-order Bessel function of the 
first kind OTT . Notice that the above resuh is identi- 
cal to that obtained in the absence of the high-frequency 
magnetic field [see Eq. (11001 )1. but with the perpendicu- 
lar frequency cl>j_ renormalized to a value (L)j^Jo(r). The 
solution of Eq. ( fT2b with the initial condition ( fTOl l is 



G(t) = e-"""'e 



-''/idr„-'fl(0)o-;/2 



(14) 



Finally, according to Eq. (|7]i, the time evolution oper- 
ator obtained by the average method is 



t/av(0 = e 



-i[t-9(t)]&J2^-ilh,„ ^-ie(0)&J2 ^ 



(15) 



3. Multiple scale method 



In this section, we will use a multiple scale method 
to obtain an approximate solution of Eq. ^ valid in a 
wider range of time scales than the one obtained by the 
averaging procedure [see Eq. (fTSlll. As will be seen in 
Sec. 3] this method provides a better understanding of 
the system dynamics than the averaging method, above 
all in the long-time limit. Since we are interested in 
the high-frequency limit Ohf » 1 , we will choose e - 
f^HF ^ 1 ^s the small dimensionless parameter in the 
multiple scale asymptotic expansions. 

In order to apply the multiple scale method, we will 
introduce the rapid time scale 



as well as the time scales 



(16) 



(17) 



with k = 1, 2, Notice that fi = f is the time 

scale in which the rotating magnetic field evolves. For 
ease of notation, a function / of all the time scales t^, 
with k > n, will be denoted by f(T„), [i.e., f(Tn) = 
f{t„,t„+u ■■■)]■ 

Even though the exact solution of Eq. ^, G{t), is 
a function of t alone, we will seek solutions of this 
equation which are functions of all the time scales in 
Eqs. (fT6b and ( fTTl l treated as independent variables. In 
other words, we will seek solutions of Eq. (|9]l of the 
form^(7"o). This is just an artifice to remove secular be- 
haviors and, at the end of the calculations, all these time 



scales will be expressed in terms of f by using Eqs. dTSI l 
and([T7]i. 

Making formal use of the chain rule for partial dif- 
ferentiation in Eq. (|9]i, it is easy to see that the operator 
Q{Tq) fulfills the differential equation 






dtk 



^ehito)giTo). 



(18) 



If one assumes a perturbation expansion of the form 

oo 

&iTo)^J]^g^"\To), (19) 



n=0 



replaces it in Eq. ( fTST l. and equates the terms with the 
same powers of e, one obtains the hierarchy of equations 



. g^<°'(To) ^ 

I 7. 0, 

dto 

for n - 0, and 



" 5^<«-*>(ro) 

k={) 



'Z 



dh 



kto)@^"-'\n) 



(20) 



(21) 



for n > 1 . This hierarchy of equations will be the start- 
ing point of the rest of our analysis. 

3.1. Time evolution operator for time scales t <K e^^ 

According to Eq. (l20l i. ^*°*(7o) does not depend on 
the rapid time scale to, i.e.. 



.(0). 



0'">(To) ^ g'.'iTi), 



^(0)/ 



(22) 



where 0\ (Ti) is an arbitrary function of all the time 
scales except to. 

Thus, setting n = 1 in Eq. (l2Tl i. one obtains 



dto 
whose solution is 

3(1) 



dti 



+ Kto)g'"iri), (23) 



&''>(To) = gy'iTi)-to 



i(i)^ 



r-to 

I dt'^kt'o 
Jo 



dti 
)Qf{Ti), (24) 



^(i)/ 



with @\ {T\) being another arbitrary function of all the 
time scales except to. 

In order to highlight the secular terms appearing in 
the above expression (i.e., the terms which diverge as 
to — » -Hoo), it is convenient to introduce the operator 



l:(fo)= I df^ [/!(?;)- ^eff] , 



(25) 



where h^ff is the effective Hamiltonian defined in 
Eq. ( [13] ). Then, Eq. ( l24l i can be rewritten in the form 



&\'\Ti) ■ 



■itito)gf\Ti) 



-to 



dgf\Ti) 



.(0)/ 



dt. 



i/!eff^r(^l) 



(26) 



From its definition, it is easy to see that £(fo) is a 
periodic function of to of period 2n and, consequently, 
remains bounded in the limit fo — > +°°- Thus, in or- 
der to avoid the secular term appearing in Eq. ( |26] |. it is 
required that 

' ^ -ih,sgf\Ti). (27) 

oti 

The solution of the above equation is 



^(0) 



^r(r,) = e-'""-^ "'(r2), 



''i'!cir/5(0)/ 



(28) 



where §2 C^i) is an arbitrary function of all the time 
scales tk with k > 2. With this choice of ^?\Ti), 



the secular term appearing in Eq. (|26] | vanishes, so that 

^<i>(ro)isgivenby 



Ml) 



(0)/ 



@''\To) = g\"iTi) - mt,)g'^'iTi) 



(29) 



Now we apply the result given in Eq. ( 128b to time 
values such that f <K e '. In this case, since tk <i: 1 for 
k>2, the function §2 vTi) appearing in Eq. ( l28T l can be 
replaced by a constant to be determined from the initial 
condition in Eq. (fTOl i. Then, making t\ = f in Eq. ( l28l l 
and taking into account Eqs. ( fT9l l and (l22T i. it results that 



G(0 = e-'^^ife 



-itK,-i0(O)&,l2 



+ 0{e). 



(30) 



This last approximate expression coincides with that 
obtained in Sec. |2]by using the averaging method. As 
we will see in the following subsections, the multiple 
scale method allows us to go beyond this result and ob- 
tain approximate expressions of the evolution operator 
valid in a wider range of time scales. 

3.2. Time evolution operator for time scales t <K e"^ 

Setting n = 2 in Eq. (|2B and using Eqs. (|22|, ^ 
and ( |29] l. one obtains 

I = - ' T(Ti ) + I U(to) Q\ '(Ti) 



dto 



+kto)0?iTi), 



where 



r(i\)- H 



dt\ 



dt2 



(31) 



(32) 



and 

n(fo) = tito)k« - kmito). (33) 

The integration of Eq. ( ISTT i results in 

^^'HTo) = #i''(7-i)+ J°df^ [tV{t'o)0f\Ti) 

-ikt[;)S^'\Ti)\- tof{T,), (34) 

with Q\ iT\) being an arbitrary function of all the time 
scales except fo- Introducing the operators 

-In 



1 r" 

Heff = — dfon(fo) 

2n Jo 



(35) 



and 



A(fo)= r°df;[n(f;)-neff], (36) 

and using Eq. ( 1251 ), Eq. ( l34l i can be rewritten in the form 

#'^(To) = gfhTi) - tof'iTx) + k{to)gf{Tx) 
-itito)&\'\Ti), (37) 



with 



+ ik«&?\Ti). 



(38) 



As in the case of the operator £(?o), from Eq. ( l36l l it 
is easy to see that A(fo) is a periodic function of fo of pe- 
riod 2n and, consequently, remains bounded in the limit 
fo — > -Hoo. Thus, in order to remove the secular term ap- 
pearing in Eq. ( l37b one has to impose the condition that 
'F'iTi) = or, taking into account Eqs. ( |28] ) and ( l32l l. 
that 



dti 



= -iks@\'\T,)-e-'"'-' 



+ n,ffe-'"'-'gf\T2). 



dgf\T2) 



dto 



(39) 



The formal solution of the above equation is 



#/'(T-i) 



gi..,r,,-„.<''^'' 



dt2 



+ k(t,)Qf{T2) 



(40) 



where 



A(f 1 ) = r ' df ; e'''i '''^ fleffe""'' '''" (4 1 ) 

Jo 



and Oy^i) is an arbitrary function of all the time 



^(1)/ 



scales except fo and t\ . By choosing Q\ {T\) of the form 
given by Eq. ( |40] |. we achieve to remove the secular term 
appearing in Eq. dJTl l. so that it reduces to 



V2) 



(0)^ 



-it(to)g\'\Ti), (42) 

In order to analyze the behavior of the operator A(f i) 
in the limit fi — > +oo, it is necessary first to evaluate the 
operator Heft- From Eqs. ( l33l l and (|35] |, together with 
Eqs. (fTTI) . ( fT3] l and dZSl ). we obtain after some lengthy 
calculations that 

-2ff 



1 r r 

fteff = -— I dfo I d?;, [/i(fo),/z(fo)] 



= Ileff ■ 0-, 



(43) 



where & - &xex + o-yCy + d-,ez and Ileff is the 3-D vector 
Ileff = T^ {(1 + i^ii) [a(r, i^)ex - b(r, i^)ey] 



■(Oj_Jo(r)a(r,(p)ei} . 



(44) 



The functions a(r, (p) and fe(r, ^) appearing in the above 
expression can be written in the form 



a{r, (f) - 2 



and 



b(r, (f) — 2 



Jo 



sin[r sin(fo)] - nHoir) (45) 



cos[rsin(fo)] - <fiJa(f) 



(46) 



with Ho(r) being the zeroth-order Struve function ||37 I. 
Once we have calculated Ileff, the operator A(ff) can 
be obtained by using the following general procedure. 
Let us introduce the effective frequency 



fieff = ^/oT^^i^^Ti^^^wowF, 



(47) 



and assume for the time being that Qeff i^ 0. Then the 
effective Hamiltonian /leff can be expressed as 



h, 



a 



■eff 



eff ■ 



no-. 



where n is the 3-D unit vector 

n = Q;^ [cu_LMr)ex + (1 + W||)ez] . 



(48) 



(49) 



If we consider now an operator of the form a • &, 
where a is an arbitrary 3-D vector, then it is a well- 
known result ll34ll that 



e"''''>^a-&e-"''''" 



aR(ri) ■ &, 



(50) 



where 

aR(ri) = (n-a)n-i-cos(nefrfi)[a-(ii-a)ii] 

H-sin(Qefffi)n A a (51) 

is the 3-D vector generated from a by a counterclock- 
wise rotation of an angle Qefff t about the unit vector n, 
and A denotes the cross product of two 3-D vectors. Ac- 
cording to Eqs. ( l50l ) and ( fSTI ). we obtain that 



Jo 



' ■^"''''"'a ■ o-e"'''i''" 



df'i e' 



af(fi)-o-. 



(52) 



with 



/. ^ / ^ . Sin(nefffl) , , , -, 

af(ff) - (n ■ a)n fi h — [a - (n ■ a)n] 



^eif 
1 -COS(Qeff fl) 



a 



n A a. 



(53) 



■eff 



Applying the above results to the case a - Ileff, and 
taking into account that the 3-D vectors n and Ileff are 
perpendicular, we obtain that 



A(fi) 



n 



eff 



sin(Oefffi) 

^eif 
1 -COS(Qeff ff) 
f^eff 



n A Ileff 



(54) 



From Eq. ( |54] | it is clear that, for Qeff ^ 0, A(fi) is a 
periodic function of f i of period Qeff and, consequently, 
remains bounded in the limit ff — > +co. As just men- 
tioned, this is a consequence of the fact that the 3-D vec- 
tors n and Ilefr are perpendicular. In the case Qeff = 0, 
[i.e., if W|| = -1 and Jo{r) - 1.38.1 1. both operators 
/jeff and rieff vanish identically and, therefore, so does 
A(fi). Thus, for any value of Qeff the operator A(fi) re- 
mains bounded in the limit fi — > H-oo and, in order to 
remove the secular term appearing in Eq. ( l40l i. one has 
to impose the condition that 



dgf\T2) 



0, 



dt2 
or, equivalently, that 



@f\T2) 



: #f (T-3), 



(55) 



(56) 



where ^3 (7^3) is an arbitrary function of all the time 

scales except to, t\ and t2- With this choice of ^2 iTi), 
the secular term appearing in Eq. ( |40] | vanishes, so that 
#/'(ri)isgivenby 



#/'(ri) = e-'"''- [g^^\T2) + A(fi)^f (71)1 ■ 



(57) 



Replacing Eq. (|56] | in Eq. (l28l l and using Eq. (l22l i. we 
obtain that 

^'"'(T-o) = Gf\Ti) = e-'"^-^f (Tl). (58) 

From the above expression and by a similar reasoning as 
the one appearing below Eq. ( |29] |. one can see that the 
approximate solution (|30] |. which was obtained under 
the assumption that f <k e ', remains valid in the wider 
range of time scales defined by the condition t <s: e^^. 

3.3. Time evolution operator for time scales t <K e"^ 

Setting n = 3 in Eq. (EJ and using Eqs. (|29]i, (HTt . 
(I42l i, ( l57T i and (|58]), one obtains that 



dto 



+ iU(to)g\'\Ti)-r(to)gf\Ti), 



where 



3 a/^(3-*)/rr \ 



jt=i 



5fi 



(59) 



(60) 



and 



r(fo) = A(fo)/zeff - h{to)A(to) + 2;(fo)neff. (61) 

The integration of Eq. (|59] | yields 



.(3)/ 



^^''(T-,,) = ^i^'(ri) 



-/ J" df,', [/!«,) 



W?\Ti) 



+ iUit'o)@\'\Ti) - tit'o)gf\Ti)] 

-to^(Ti), (62) 

where g\ {T\) is an arbitrary function of all the time 
scales except fo- 

Taking into account that the operators h{tQ), A(fo) and 
£(fo) are all periodic functions of fo of period 2n, so will 
be the operator r(fo). Thus, if we define the operators 

reff = — I dfof(fo) (63) 



and 



-f 

27r Jo 

Y(fo) = J°df;[f(f;)-feff], 



(64) 



we can assure that T(fo) is also a periodic function of fo 
of period 2n. Replacing the definitions dZSl l, (|36] | and 
(l64l l in Eq. (l62l i. one obtains that 

^'"(T-,)) = #i'^(ri) - fo^'(T'i) - it{tQ)gf(T,) 



+ A(fo)^"(ri) + /T(fo)^"'(ri), 



(65) 



with 

^'(T-i) = ^(T-i) + //ieff#''(T-i) 

-fteff ^I'^rO - / feff^f (Ti). (66) 

Since the operators £(fo), A(fo) and T(fo) in Eq. ( l65T l 
are periodic functions of fo of period 2n, they remain 
bounded in the limit fo — > +°°- Therefore, in order to 
remove the secular term appearing in Eq. ( l65T l one has to 
impose the condition that ^'{T\) - 0. Taking into ac- 
count Eqs. (EB, (|57]l, (|58]l, (l60]l and (|66]l, the condition 
y['(7~i) = leads us to the diff'erential equation 



5fi 



£'(T2) 



-^^<"(r2)-s(f0^f(T-3) 

Ofi 



where 

and 

S(fi) = 



5^"(T2) dgxHTi) 



dto 



gA(fi) 

5f, 



5f3 



ifi/icirf 



A(fi) + /e"'«'^reife 



„-it\Ks 



(67) 



(68) 



(69) 



The formal solution of Eq. (|67] ) is 



^f(T-i) 



,-»i/ic(r r/2(2) 



^ ^^(7-2) - fi £)(r2) 



■A(fi)^^'^(r2) + ©(fi)^f(r3)], 



(70) 

where @2 i'^i) is an arbitrary function of all the time 
scales except fo and t\, and 

0(fi)= r'df;s(f;). 

Jo 



(71) 



In order to know the behavior of the operator 0(fi) in 
the limit fi — > +oo, we first have to evaluate the operator 
Feff. Using Eq. ( |6TI ). together with Eq. ( |25] ), it is easy to 
show that 



feff = [Aeff,/leif] + ^ f dfo S(fo)fl(fo), (72) 



with 



1 r^'^ . 

Aeif = — dfo A(fo). 

27r Jo 



(73) 



Replacing in Eq. ( l72l l the definitions ( |25] | and ( |33] ). and 
using Eqs. (fTTI) and ( fT3T l. one obtains after some lengthy 
calculations that 



f eif = [Aeif , /leff] + Q ' O", 



(74) 



where 



and 



q = - (^) [y- [«-^('-' ^)ex + ay(r, ^)ey] 6(fi) = / J ' df; e"'.''- [Aes, ^eff] e-"'.''- 



+ (1 + W||)a',(r, <^)e2 



with 



ax{r,ip) = ^— [a(r,(^)]^ - yi(r), 



a,.(r, ^) 



2 
/o(r) 



fl(r, (^) Z7(r, ^) 



and 



a,(r, (p) - 



a(r, if) 

2 


2 
+ 


b(r,ip) 

2 



- nir). 



(75) 

(76) 
(77) 

(78) 



In the above expressions yi(r) and ')'2('") are two func- 
tions of the parameter r which are given by 



7i(r) 



■jQ(r)[Ho(r)f 



2 r c^ c^ 

+ — I dip I dip' I dip" sin(r sin ip) 
^ Jo Jo Jo 

xcos(rsin(^')sin(rsin^") (79) 



and 



7T An 

72{r) = ^ [HQ{r)f - — [Mr)f 



Mr) 



dv 





27r JO 



dipip cos(rsin(^) 



1 r c^ 

— I dip I d(^'(^ cos[r(sin(^ - sin(/3')]■ 
7!" Jo Jo 



(80) 

As in the previous subsection, now we will consider 
separately the cases Qeft + and Qeff = 0. For QefF + 0, 
the first term of the right-hand side of Eq. ( |69] l can be 
easily evaluated using Eq. (l54l i. The second term can be 
calculated from Eq. ( l74l l. making use of Eqs. (ISOl l and 
( BTI ) with a = q. Replacing the results so obtained in 
Eq. (ItTI) . one finally finds that 



0(fl) = -/?7fl/left+<l)(fi), 

where 

2 / Ileff ■ Ileff 

??= 75- nq 



(81) 



2 I Qeft/ 



2r 2 



r-/o(r)ri('-) + (i + wii)'r2(r) 



(82) 



-I-; 



sin(Qefffi) 

Qeff 
1 - COS(Qeff fl) 
Qeff 

-I- / (n A q) ■ (T 



[q ■ d- - 77 /zeff] 



n q 



//Qeff 



(83) 



Notice that the first term of the right-hand side of 
Eq. ( l83T l is a periodic function of fi of period Qeff. In 
order to see that this is so, one just has to apply the same 
reasoning as for the operator A(f 1 ), taking into account 
that, as can be easily shown, [Agff , /leff ] - / p ■ o" with p 
being a certain 3D-vector perpendicular to n. Accord- 
ing to Eq. ( l83T l. we can therefore assert that <l>(fi) is a 
periodic function of t\ of period Qeff. 
Replacing Eq. ( ISTT i in Eq. ( iTOl i. it results 

gf{T{) = e-''>''"^ \9f{Ti) - h i)'{T2) 

+ k{t,)g^^\T2) + Hti)Qf{T3)\, 



where 



i)'(T2) = D{T2) + iiik«Gf{Ti). 



(84) 



(85) 



Since both A(/i) and <l)(fi ) are periodic functions of t\ of 
period Qeff, they remain bounded in the limit fi — > +00. 
Therefore, in order to remove the secular term appear- 
ing in Eq. ( l84b one has to impose the condition that 
D'{T2) - or, equivalendy, that 



dt2 



' +ir]h,sgf\T,) 



dh 



(86) 



The integration of the above equation leads to 



G^2\T2) = S^i\T,) 



■h 



' +//7/leff^f(T3) 



dh 



(87) 



where Q^^ (T3) is an arbitrary function of all the time 
scales except fo, fi and f2- 

To avoid the secular term appearing in Eq. ( l87l l we 
have to assume that 






= -iT]kff&f\T,), 



(88) 



whose formal solution is 



^(0) 



iniihcir AW ( 



GV^T-i) = e-""-'"-^7(r4), 



(89) 



^(0) 



with Q\ {Ti) being an arbitrary function of all the time 
scales except to, t\, ti and f3. Replacing the above result 
in Eq. ^, we finally obtain that, for Qeff ^ 0, ^*"'(To) 
is given by 



g(^\To) = e-'("^''«''-^f (r4). 



(90) 



Let us consider now the case in which r is a zero of 
the Bessel function Jo{r) and (jl)\\ - -1, so that Qeff = 0. 
Specifically, we will assume that r - rj, where rj is the 

y'-th zero of Jo{r) with y = 1,2, As mentioned in 

the previous subsection, in this case the operators /letf, 
rieff and A(fi) vanish identically. Thus, from Eqs. (|69] l, 
(TTTT i. ( I74l i, ( |75] ). ( l76b and ^7} it is quite easy to see that 



0(fl) = /feff fl = M -T^j TlJ^I O'x' 



(91) 



where yiy = yiifj), i.e.. 



^ Jo 



I di^' I d(^" sin(rj sin ^) 
Jo Jo 

X cos(rj sin ip') sm(rj sin tp"). (92) 

Replacing Eq. ( |9T1 ) in Eq. ( TTOb . and taking into ac- 
count that A(f 1 ) = 0, one obtains that 

gf\Ti) = e-"''''" [gf(T2) - h D'(T2)] , (93) 
with 

D'(T2) = D(T2) - i (^) 71.J &, gfiTi) . (94) 

Then, following the same procedure as in the case Qeff i^ 
0, it results that 



0^°\To) = ^^(73) = e'(^y^^-^^''^gf\T4) 



(95) 



when Qeff = 0. 

Now we apply the same reasoning as at the end of 
subsection 13.11 Let us assume that we are only inter- 
ested in time values such that f «: e^^, so that, according 
to Eq. ([nil, it is fulfilled that f^ <k 1 for A: > 4. Then, in 
Eqs. ( l90b and (l95T l, it is possible to replace the function 
^4 (T4) by a constant to be determined from the initial 
condition in Eq. ( fTOt . Expressing f 1 and t^ in terms of t 
and taking into account Eq. ( fT9l ). we finally obtain that 

G(t) = e-''''n„g-'e(0).^./2 + ^(^)^ (96) 

where h^^ is a multiple scale Hamiltonian defined by 



(1 +e^ri)h 



eS 



(97) 



in the case that Qeff ?t 0, and by 

v3 



h„ 






(98) 



in the case that Qeff = and r = rj. 

Finally, according to Eq. d?), the time evolution oper- 
ator obtained by the multiple scale method and vaUd for 
time scales t «; e^^ is 

U^M = e-'['-fiW]<^./2^-,Y/,„,,^-,-fl(0R,/2^ (99) 

The above three expressions, together with the defini- 
tions of 77 [Eq. ( l82] ll and jij [Eq. (l92] il are the main 
results of this section. 



4. Results 

We begin this section with a brief description of the 
resonant behavior displayed by a spin- 1/2 particle under 
the influence of the rotating magnetic field in Eq. ([TJ, 
when the high-frequency magnetic field in Eq. (O is ab- 
sent (i.e., when r = 0). Then, according to Eqs. dHJ and 
(fTTT i. the Hamiltonian /^(QhfO becomes time indepen- 
dent and equal to 



^0 = -J [f^J-O-x + (1 + W||)o-.] . 



(100) 



Taking into account Eqs. (|7]i and (|9]l, the time evolution 
operator will be then given by 



U(t) = e 



-kd-;/2 -itho 



(101) 



Henceforth, we will be interested in studying the time 
dependence of the expected value of &,. If we take as 
initial condition an arbitrary state \^, 0), this expected 
value can easily be evaluated from the time evolution 
operator by using the expression 



{&,;t} = {'¥,0\U''{t)&M{t)\^',Q). 



(102) 



Let us assume that we take as initial condition I^F, 0) - 
|±), with |±) being the eigenvectors of cr- corresponding 



to the eigenvalues ±1. Then, from Eqs. ( llOlb and ( I102l i. 
together with Eqs. ( fSOl l and ( fSTT l with /leff replaced by /jq 
and f 1 by t, it is easy to see that 



w. 



(a-,; t^^^ ^±\l + ^ [cos(Qor) - 1] 



(103) 



where <o-,;f)(±> = {±\U\t)&,U(t)\+} and Qq = 



^ 



(1 -H W||)2 + aj]_. 

From the exact expression ( 1103b one can see that, 
in the absence of the high-frequency magnetic field, 
(o-^; f)'=^' is a periodic function of time of period 27r/Q(). 



The average of this function over a period is ±(1 + 
cj\\)'^IQ^y and the ampUtude of the oscillations is given 
by 






(1 + (y||)2 + ijj\ 



(104) 



This amplitude as a function of co\\ displays a resonant 
behavior with a maximum at wy = - 1 . This behavior is 
depicted with a solid line in Fig.|2]for the case cj^ - 3. 
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Figure 2: Resonant behavior of the amplitude A(w||) for Wx = 3 and 
r = 0, 1 and 2. Solid line: exact result for the case r = obtained 
from Eq. 1104) . Dashed and dotted lines: analytical results obtained 
from Eq. 11081 for r = 1 (dashed line) and r = 2 (dotted line). Sym- 
bols: results obtained from the numerical solution of Eqs. (4)-(6), and 
the use of Eq. 1102) with I'F, 0) = |+). The parameter values of the 
high-frequency magnetic field used in the numerical calculations are: 
^HF = 50, If = nil, r = 1 (crosses) and r = 2 (pluses). 

When the high-frequency magnetic field is present, 
we do not have an exact expression for the time evo- 
lution operator and have to resort to Eqs. (fTsT i or (|99] |, 
depending on whether one uses the averaging method or 
the multiple scale method. Let us consider first the case 
in which r is not a zero of J(,{r), so that Qgff + 0. Then, 
from Eqs. ([T3ll, ([TSll and ( 1102b . together with Eqs. (ISOl l 
and (ISTl i. the averaging method yields 



<<^z;o: 



(±) 



± 1 + 






[cOS(QeffO - 1] 



(105) 



However, if one uses the multiple scale method, from 
Eqs. ( |97] ) and ( |99] ) one obtains 



/ ' a(±) ^ 1 1 , [^-L-^o('-)] r /o 

{o-z; trj = ± n + — -5 — [cos(n„ 

I ^ eff 



.0-1] 



with the multiple scale frequency 

Qms = (1 +e^7;)Qeff. 



(107) 



Notice that both methods lead to the same expression 



[<o±Jo(r)f 

"eff 



[c^±Jo(r)Y 



(i + wii)2 + [wx/o('-)r 



(108) 



for the amplitude of the oscillations of the function 

(o-,;f)<±'. 

-1 



(106) 



According to Eq. dlOSb . the resonant peak at wy 
persists in the presence of a high-frequency magnetic 
field when r is not a zero of Jo(r). Since the width of 
the peak is related to co^Joir), the presence of the high- 
frequency magnetic field gives rise to a narrowing of 
the resonant peak with respect to the case r - Q. These 
analytical results are shown in Fig. |2] for cl>± = 3 with 
a dashed line (r = 1) and a dotted line (r - 2). In 
Fig. 121 we also show the amplitude values obtained from 
the numerical solution of Eqs. (|4li-(|6]), and the use of 
Eq. (11021 ) with I^F, 0) = \+). The numerical integration 
of the differential equation (|4| has been performed using 
the function NDSolve of the software package Mathe- 
MATicA with the option MaxSteps — > 00. In order to have 
a well-defined amplitude value, the high-frequency and 
small-amplitude oscillations observed in the numerical 
values of {&,; f)^^' have been removed by performing a 
time average of {&~; f)*^* over a high-frequency period. 
The parameter values of the high-frequency magnetic 
field used in the numerical calculations are: Q.uf - 50, 
ip = 7t/2, r - 1 (crosses) and r = 2 (pluses). Notice that 
there is a very good agreement between the analytical 
and the numerical results. 

Although the analytical expressions for {&,; f)*-"^' ob- 
tained by the averaging method [Eq. ( 11051 )1 and the 
multiple scale method [Eq. ( 11061 )1 look very similar, 
there exists a slight difference Qm, - Qeff = e^yyQeff 
between the frequencies of the oscillations predicted 
by both methods. This slight difference becomes im- 
portant for large values of t. In particular, for time 
scales t x n/\£ rjQ.ef{\, the expected values (o-,; f)m,s ^nd 
((T,; f)av ^re n radians out of phase with each other. For 
instance, for a>± - 3, (i)\\ = -I, r = I and Qhf - 50, this 
situation occurs for t x 7r/\£^riQ.gff\ » 1700. To work out 
thisnumber, we have used the value yi(l) ^ -0.684533, 
obtained numerically from Eq. ( |79] l, for the calculation 
of T] [see Eq. (l82bl. The integral appearing in the defi- 
nition of yi(l) [see Eq. ( l79b l has been evaluated numer- 
ically using the function NIntegrate of the computer 
algebra program Mathematica. 

In Fig. [3l it is shown the time dependence of 
(d",;?)"-^' in the intervals t € [0,6] (upper panel) and 



10 



t e [1700, 1706] (lower panel), for the above men- 
tioned parameter values. The results obtained from 
Eq. ( 11051 ) (averaging method) have been plotted with 
dashed lines, and those obtained from Eq. ( 1106b (multi- 
ple scale method) with solid lines. We also depict with 
crosses the results obtained from the numerical solution 
of Eqs. &-© (for (p = 7t/2), and the use of Eq. ( fT02l l 
with I*?, 0) = |-i-). For short time scales (upper panel), 
both the averaging and the multiple scale methods agree 
very well with the numerical results. In fact, the dashed 
line is completely hidden by the solid one. Neverthe- 
less, in the time interval [1700, 1706] only the multiple 
scale method reproduces the numerical results, whereas 
the averaging method predicts a time dependence which 
is approximately n radians out of phase with respect to 
these ones. 
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Figure 3: Time dependence of (&-;t)^^'' in the time intervals t e 
[0,6] (upper panel) and t e [1700,1706] (lower panel). Dashed 
lines: analytical results obtained from Eq. 11051 (averaging method). 
Solid lines: analytical results obtained from Eq. )106t (multiple scale 
method). Crosses: results obtained from the numeiical solution of 
Eqs. E)-|6j, and the use of Eq. fT02t with I'P, 0) = |+). The parameter 
values are: cj^ = 3, mn = -1, r = 1, Ohf = 50 and tp = n/2. In the 
upper panel the dashed Hne is hidden by the solid line. 

Now let us consider the case in which r is equal to one 
of the zeros of the Bessel function Joir), for instance, 
the j-th zero rj. If we also assume that W|| i^ -1, it is 
easy to see that Eqs. dlOSb and (1106b remain valid, since 
the condition that Qetj i^ still holds. Consequently, 
from Eqs. (1105b and ( 1106b one obtains that 

{&,;ty£>^{&,;tt!^±l, (109) 

for r - rj and wy it -1. 

For r - rj and wy = -1 the averaging and the multiple 
scale methods lead to different expressions for {&~; f)'=^\ 
If one uses the averaging method and takes into account 
Eqs. ( fTSl l and (1102b . the result is (o"c;;f)av = ±1^ since 



the effective Hamiltonian in Eq. ( fTlT l vanishes when 
r - rj and wy - -1. By contrast, if one uses the mul- 
tiple scale method and takes into account that for these 
parameter values fleff = 0, from Eqs. ( l98T l and (l99T l one 
obtains 



i^z-J) 



(±) 



±cos 



(Qms,/) , 



with the multiple scale frequency 



Q 



msj 



j'^^rLj- 



From these results one concludes that, for r 



(110) 



(111) 



the 



averaging method predicts that Aav(w||) = for all the 
values of wn and, consequently, the complete disappear- 
ance of the resonant behavior By contrast, the multiple 
scale method predicts that 



An,s(a)\\) = 



if co\\ + 
if iD\\ - 



(112) 



i.e., the resonant peak at wy = -1 persists but its width 
becomes zero. This last behavior is depicted in Fig. |4] 
with a dashed line. In Fig.lH we also show with crosses 
the amplitude values obtained numerically, by the nu- 
merical procedure described previously, for the param- 
eter values (1)^ - li, r - r\ ■x, 2.40483, Qhf - 50 and 
ip - 71 12. The numerical results highlight the existence 
of an extremely narrow peak around wy = - 1 and, con- 
sequently, confirm the prediction made by the multiple 
scale method. Nevertheless, the multiple scale method 
is not able to explain the finite width of the resonant 
peak observed in the numerical results shown in Fig. |4l 
This fact shows up the inability of the analytical expres- 
sion in Eq. ( 11091 ) to describe properly the time behavior 
of ((T;;f)*''* in a small neighborhood around wy = -1. 
This inability is not surprising if one takes into account 
that the expression for (o-j,;?)!,^ in Eq. (11091 ) does not 
converge to the result in Eq. ( II 10b in the limit co\\ — > -1. 
In order to avoid this discontinuity, we suspect that it is 
necessary to extend the multiple scale method described 
in Sec. [3] to a larger number of time scales. 

To illustrate this point, in Fig.|5]it is shown the time 
dependence of {&/, f)*^* for r - ri ^ 2.40483 and five 
values of wy around u)\\ - -1. The rest of the parame- 
ter values are: u)_i_ - 3, Qhf = 50 and ip - jijl. The 
analytical results obtained using the averaging method 
[i.e., Eq. ( 11091 ) for all the values of wy] are plotted with 
dashed lines, and those obtained from the multiple scale 
method [i.e., Eq. ( |109b for wy i^ -1, and Eq. ( II 10b 
for wy = -1] with solid lines. We also show with 
symbols the results obtained numerically. In the up- 
per panel, it is considered the case wy = -1.1. In this 
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Figure 4: Resonant behavior of the amphtude A(w||) for r = ri x 
2.40483 [first zero of the Bessel function Joir)]. Dashed line: analyti- 
cal results obtained by the multiple scale method [Eq. )112H . Crosses: 
numerical results obtained following the same procedure as in Fig. [2] 
The parameter values used in the numerical calculations are: (H^ = 3, 
f^HF = 50 and tf = n/2. The scale of abscissa has been magnified 
with respect to Fig. |2]in order to make visible the finite width of the 
resonant behavior observed numerically. 



case, the value of wy is not too close to the resonant 
value wy - -I, and both the averaging and the multiple 
scale methods reproduce properly the numerical results 
depicted by crosses. In the middle panel, it is consid- 
ered the cases wy = -1.005 (pluses), wn = -1.001 (cir- 
cles) and a>\\ = -1.0005 (triangles). As a>\i gets closer 
to the resonant value wy - -1, the disagreement be- 
tween the analytical results in Eq. 



I109l l and the numer- 
ical ones becomes more pronounced, and the numeri- 
cal results converge gradually to the multiple scale so- 
lution in Eq. dllOb . Since the averaging and the multiple 
scale solutions coincide when wy ^ -1 [see Eq. ( I109H , 
the dashed lines are hidden by the solid lines in the up- 
per and middle panels. Finally, in the lower panel it is 
shown the case wy = -1. In this case, only the multiple 
scale expression in Eq. 



I ll 10b is able to reproduce the 
numerical results depicted by squares. To calculate the 
multiple scale frequency in Eq. ( Illlb . we have used the 
value 7ij « -0.603984, which has been obtained nu- 
merically from Eq. ( |92] | using the function Nlntegrate 
of the program Mathematica. 

We finish this section with a brief discussion of the 
role played by the phase shift ip between BHF(f) and B(f) 
in the system dynamics. Since the frequencies Qgff , f^ms 
and Qms,/ do not depend on ip, the analytical results ob- 
tained up to now in this section are also independent of 
this parameter. As we shall see below, this is so due to 
the initial conditions considered, namely |*P, 0) = |±). In 
the numerical results, the value of ip affects only slightly 
the high-frequency and small-amplitude oscillations ob- 
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Figure 5: Time dependence of (iT;;f)'^' for r = r[ x: 2.40483 [first 
zero of the Bessel function Joir)] and five values of wn around wy = 
-1. The rest of the parameter values are: 6Jx = 3, flnp = 50 and 
tfi = n/2. Dashed lines: analytical results obtained from Eq. U09t for 
all the values of wy (averaging method). Solid lines: analytical results 
obtained from Eq. 1109) for to^ # -1, and Eq. jllOt for ajy = -1 
(multiple scale method). Symbols: numerical results. Upper panel: 
<jj|l = -1.1 (crosses). Middle panel: wy = -1.005 (pluses), niy = 
-1.001 (circles) and(ij|| = -1.0005 (triangles). Lower panel: wy = -1 
(squares). In the upper and middle panels the dashed lines are hidden 
by the solid lines. 



served in {o-.; t}^^^ (see Fig. |6]l. 

A more interesting situation arises when one consid- 
ers an arbitrary initial condition of the form 



l^V, 0) = c 1+) + e'"* Vl-c2|-), 



(113) 



where c e [0,1] and e [0, 27r). In this case, it is 
straightforward to see from Eqs. ( |99] | and ( 1102b that, if 
we write explicitly the dependence of the expected value 
of (Tj on the parameters if and (p, it results that 



(d-z;t,{<f, 



{&,-j,{0,4>'])„ 



(114) 



where 0' - (p + rsin (f. That is, in order to calculate 
the expected value of o-j. by the multiple scale method, a 
high-frequency magnetic field with a phase shift (p can 
be replaced by another one with a phase shift equal to 0, 
if the phase of initial state is also replaced by another 
value (p' - (f) + rsimp. From Eqs. ( fTSl ) and ( 11021 ). it 
is easy to see that this result also holds if one uses the 
averaging method. 

As an example of this property, let us consider the 
case in which r = rj and wy = - 1 . For these parameter 
values, the effective Hamiltonian in Eq. ( fT3l l vanishes 
and, consequently, from Eqs. (fTSl ), ( 11021 ) and ( 11131 ). the 
averaging method yields 



<o-,; f)a 



2c^-l, 



(115) 
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t 

Figure 6: Time dependence of {&-;!)'-•*'' for (tiy = -1, aix = 3, ilnF = 
50, r = ri X 2.40483 [first zero of tlie Bessel function /o('')]> 9 = 
njl (upper panel) and ^ = (lower panel). Dashed line: analytical 
results obtained by the multiple scale method [Eq. UIOH . Solid lines: 
numerical results. 




4000 



Figure 7: Time dependence of {&-; t) for aij| = -1, Wx = 3, Ohf = 
50, r = ri X 2.40483 [first zero of the Bessel function /o('')] and 
an initial condition of the form 1113) with c = 2"''^. Upper panel: 
</; = ff/2 and = 0. Lower panel: ^ = and <f> = ri. Dashed 
lines: analytical results obtained from Eq. U15t (averaging method). 
Solid lines: analytical results obtained from Eq. )116t (multiple scale 
method). Crosses: numerical results. 



which is independent of ip and (p. By contrast, if one 
uses the multiple scale method, from Eqs. (|98] l, (|99] l, 
(IIO2I 1 and ( III3I 1, one obtains 



(o-,;f)ms = (2C^- 1)C0s(Q™s,/) 



■ 2c Vl -c^ sin (0') sin Ul^^^A , 



(116) 



which depends on ip and 4> through 4>' ■ 

In order to compare these analytical expressions with 
the numerical results, in Fig. Qit is shown the time de- 
pendence of {&-; t) for the following parameter values: 
W|| = -1, wx = 3, Ohf = 50 and r = ri » 2.40483. In 
the upper panel, it is considered a phase shift ip - njl 
and an initial condition of the form ( 11131 ) with c - 2 '^^ 
and (p - 0. In the lower panel, it is assumed a phase 
shift if - Q and the same initial condition as before 
but with (p - r\. According to Eq. ( 1116b . the multi- 
ple scale method leads to the same conclusion in both 
panels (solid lines). The analytical results obtained us- 
ing the averaging method [Eq. (I115b l are plotted with 
dashed lines. We also show with crosses the results ob- 
tained numerically. Once again, only the multiple scale 
method is able to reproduce the numerical results, with 
exception of the high-frequency and small-amplitude 
oscillations. 

5. Conclusions 

We have studied how the resonant behavior displayed 
by a spin- 1/2 particle under the influence of a rotating 
magnetic field is aff'ected by the application of an addi- 
tional high-frequency magnetic field along the rotation 



axis. Analytical expressions for the time evolution op- 
erator and the expected value of &, have been obtained 
by using two alternative methods, namely, the averaging 
method and the multiple scale method. 

When the ratio of the strength of the high-frequency 
magnetic field to its frequency (in dimensionless units) 
is not a zero of the Bessel function J{)(y), both meth- 
ods lead to the same conclusion about the effect of the 
high-frequency magnetic field on the resonant behavior. 
Specifically, its presence only gives rise to a narrowing 
of the resonant peak. By contrast, when this ratio coin- 
cides with a zero of 7o(r), both methods lead to diff^erent 
predictions. More precisely, the averaging method pre- 
dicts the complete disappearance of the resonant behav- 
ior, whereas according to the multiple scale method the 
resonant peak persists but its width becomes zero. This 
discrepancy in the conclusions of both methods is simi- 
lar to the one appearing in the study of the phenomenon 
of coherent destruction of tunneling. There, the aver- 
aging method predicts the complete destruction of tun- 
neling, whereas other more accurate methods, among 
them the multiple scale method, conclude that this de- 
struction is only temporary. In our case, the numerical 
results show up the existence of an extremely narrow 
resonant peak and, consequently, confirm the prediction 
made by the multiple scale method. An open problem 
which arises in this context is how to extend the multiple 
scale method proposed in this paper in order to explain 
the finite width of the resonant peak observed in the nu- 
merical solution. 

Comparison of the analytical results for the expected 
value of &, with the numerical ones shows that the mul- 
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tiple scale method also allows a better description of the 
time evolution than the averaging method, above all in 
the long-time limit. 

Finally, we have studied the role played by the phase 
shift between the high-frequency and the rotating mag- 
netic fields in the system dynamics. From this study it 
follows that a modification of this phase shift is basi- 
cally equivalent to a change in the phase of the initial 
state. 

6. Acknowledgements 

This paper is dedicated to Professor Peter Hanggi on 
the occasion of his 60th birthday. The author acknowl- 
edges the support of the Ministerio de Ciencia e Inno- 
vacion of Spain (FIS2008-02873) and the Junta de An- 
dalucia. 

References 



[31] 
[32] 
[33] 
[34] 

[35] 



[36] 
[37] 



[38] 



M. Frasca, Phys. Rev. B 71 (2005) 073301. 
Q. Xie, J. Phys. B: At. Mol. Opt. Phys. 42 (2009) 105501. 
J. Shao and P Hanggi, Phys. Rev. A 56 (1997) R4397. 
A. Galindo and P. Pascual, Quantum Mechanics I and II, 
Springer- Verlag, Berlin Heidelberg, 1991. 
C. M. Bender and S. A. Orszag, Advanced Mathematical Meth- 
ods for Scientists and Engineers, Mac-Graw Hill, Inc., New 
York, 1978. 

M. lanowicz, Phys. Rep. 375 (2003) 327. 
M. Abramowitz and I. A. Stegun (eds.). Handbook of Mathe- 
matical Functions: with Formulas, Graphs, and Mathematical 
Tables, Dover, New York, 1972. 

We assume that cj^ # 0, because otherwise the solution of 
Eq. (4j is trivial. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



[8] 

[9] 
[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 
[18] 

[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 

[28] 
[29] 
[30] 



I. Schwinger, Phys. Rev. 51 (1937) 648. 

1. 1. Rabi, Phys. Rev. 51 (1937) 652. 

F. Bloch and A. Siegert, Phys. Rev. 57 (1940) 522. 

H. C. Torrey. Phys. Rev. 59 (1941) 293. 

S. H. Autler and C.H. Townes, Phys. Rev. 100 (1955) 703. 

M. Grifoni and P Hanggi, Phys. Rep. 304 (1998) 229. 

T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, G. Schon and 

W. Zwerger, Quantum Transport and Dissipation, Wiley- VCH, 

Weinheim, 1998. 

F. Grossmann, T. Dittrich, P. Jung and P. Hanggi, Phys. Rev. 
Lett. 67 (1991) 516. 

R. Bavli and H. Metiu, Phys. Rev. Lett. 69, (1992) 1986. 

R. BavU and H. Metiu, Phys. Rev. A. 47, (1993) 3299. 

M. Morillo and R. I. Cukier, J. Chem. Phys. 98, (1993) 4548. 

Y. Dakhnovskii, V. Lubchenko and R. D. Coalson, Phys. Rev. 

Lett. 77, (1996) 2917. 

M. Grifoni, L. Hartmann and P. Hanggi, Chem. Phys. 217, 

(1997) 167. 

G. Delia Valle, M. Omigotti, E. Cianci, V. Foglietti, P. Laporta 
and S. Longhi, Phys. Rev. Lett. 98 (2007) 263601. 

H. Lignier, C. Sias, D. Ciampini, Y. Singh, A. Zenesini, O. 
Morsch and E. Arimondo, Phys. Rev. Lett. 99 (2007) 220403. 

E. Kierig, U. Schnorrberger, A. Schietinger, J. Tomkovic and M. 
K. Oberthaler, Phys. Rev. Lett. 100 (2008) 190405. 

F. Grossmann and P. Hanggi, Europhys. Lett. 18 (1992) 571. 

J. M. Gomez Llorente and J. Plata, Phys. Rev. A 45 (1992) 

R6958. 

L. Wang and J. Shao, Phys. Rev. A 49 (1994) R637. 

H. Cheng, M. L. Ge and J. Shao, Phys. Lett. A 191 (1994) 1. 

Y Kayanuma, Phys. Rev. A 50 (1994) 843. 
X. G. Zhao, Phys. Lett. A 193 (1994) 5. 

C. E. Creffield, Phys. Rev. B 67 (2003) 165301. 

C. E. Creffield, Europhys. Lett. 66 (2004) 631. 

D. H. Dunlap and V. M. Kenkre, Phys. Rev. B 34 (1986) 3625. 

Y Kayanuma and K. Saito, Phys. Rev. A 77 (2008) OIOIOI(R). 
J. C. A. Barata and W. F. Wreszinski. Phys. Rev. Lett. 84 (2000) 
2112. 

J. C. A. Barata and D. A. Cortez, Phys. Lett. A 301 (2002) 350. 
J. C. A. Barata and D. A. Cortez, J. Math. Phys. 44 (2003) 1937. 
M. Frasca, Phys. Rev. B 68 (2003) 165315. 



14 



